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Abstract 

In this paper we propose a new method for studying spectral properties of the 
non-hermitian random matrix ensembles. Alike complex Green's function en- 
codes, via discontinuities, the real spectrum of the hermitian ensembles, the 
proposed here quaternion extension of the Green's function leads directly to 
complex spectrum in case of non-hermitian ensembles and encodes addition- 
ally some spectral properties of the eigenvectors. The standard two-by-two 
matrix representation of the quaternions leads to generalization of so-called 
matrix- valued resolvent, proposed recently in the context of diagrammatic 
methods [H 121 EJ IH EJ E] ■ We argue that quaternion Green's function obeys 
Free Variables Calculus |Z1IH|- In particular, the quaternion functional in- 
verse of the matrix Green's function, called after [Oj Blue's function obeys 
simple addition law, as observed some time ago [HE!- Using this law we derive 
new, general, algorithmic and efficient method to find the non-holomorphic 
Green's function for all non-hermitian ensembles of the form H + iH , where 
ensembles H and H are independent (free in the sense of Voiculescu |Zj ) her- 
mitian ensembles from arbitrary measure. We demonstrate the power of the 
method by a straightforward rederivation of spectral properties for several 
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examples of non-hermitian random matrix models. 
PACS: 05.40.+j; 05.45+b; 05.70.Fh;11.15.Pg 
Keyword: Non-hermitian random matrix models. 



1 Introduction 

Random matrix models provide a powerful framework for modeling numerous 
physical phenomena, with applications covering all branches of theoretical 
physics |TTJ EH El Among different classes of random matrix models the 
non-hermitian random ensembles form a fascinating class. Contrary to the 
hermitian ensembles, where real eigenvalues form cuts on real axis, general 
non-hermitian ensembles develop spectrum, which covers two-dimensional, 
often multiple-connected support on the whole complex plane. From mathe- 
matical point of view, non-hermitian ensembles are challenging, since several 
standard methods of hermitian random matrix calculus fail in this case. On 
the other side, the issue of non-hermitian random matrix ensembles is far 
from being academic. Non-hermitian ensembles are omnipotent in various 
branches of physics and in interdisciplinary sciences. Sample applications 
are: open chaotic scattering ^3] , spectral properties of Euclidean Dirac op- 
erators in the presence of chemical potential or CP-violating angle 9 in 
Quantum Chromodynamics ^3], non-hermitian generalizations of Anderson 
localization in mesoscopic systems ^Hj, modeling of the chemical transitions 
in dissipative systems matrix generalizations of multiplicative diffusion 
processes |TH| or evolution of the spectral curves for non-hermitian ensembles 
in the context of the growth problem [T^]. All these applications call for new 
calculational methods to deal with the problem of complex spectra. 

In this paper we present a new technique, introducing the quaternion gen- 
eralization of the Green's function, the quaternion generalization of its func- 
tional inverse (quaternion Blue's function) and we formulate the quaternion 
addition law for non-hermitian ensembles. In section 2, we briefly remind 
basic known facts on hermitian and non-hermitian Green's functions. In sec- 
tion 3, we make a connection to Free Random Variables calculus. We com- 
ment also on some previous approaches based on diagrammatic techniques. 
In section 4, after establishing the notation, we outline the construction and 
we present the quaternion generalization of the hermitian Green's function 
and its properties. In particular, we present a general form of the addition 
law for non-hermitian ensembles. 

In section 5, we formulate the operational form of the addition algorithm. 
In section 6, we adapt the operational form of the addition formalism for 
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the particular case of the Gaussian randomness. Then, to demonstrate the 
power of the method, we provide straightforward derivations of three classical 
results in non-hermitian random matrix models. Section 7 summarizes our 
conclusions. Appendices hide some necessary technical details and/or proofs. 



2 Hermitian and non-hermitian Green's func- 
tions 

Let us first focus on a hermitian random matrix model. A basic tool to 
investigate eigenvalues' distribution for if is a resolvent, or Green's function, 

where ljy is a unit matrix of size N. The usefulness of this function stems 
from the fact, that, for any finite N, Gh{z) can be written as 

1 - 1 

G *W = n(Ej=x)> A ' gR ' ( 2 ) 

i=l 1 

after diagonalizing if by a unitary similarity transformation, so that Gh{z) 
is a meromorphic function with N poles at A» on the real line. In the large 
N limit the poles merge into cuts on the real axis. 

Green's function Gjj(z) can be used to reconstruct the spectral density 
function Ph(X) due to the relation 

— ^" = PV^ - m6(X) (3) 
A + le A 

where PV denotes the principal value distribution and e is meant implicitly 
approaching limit e — > 0, so that 

PhW = ^(Tr6(Xl N - H)) = --lmG H (X + ie). (4) 

Hence Ph{X) can be read out from discontinuities of the imaginary part of 
the Green's function. 

There are several ways of calculating Green's functions for hermitian ran- 
dom matrix models piUITH IT^]. Here we mention the diagrammatic approach, 
after [20] • A starting point of the approach is the expression allowing for the 
reconstruction of the Green's function from all the moments. 

n>0 
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The resolvent Gh{z) is interpreted as a series in the guise of the Feynman-like 
diagrammatic expansion in the large number of colors limit. This series can 
be efficiently evaluated exploiting the analogy of RMT to + dimensional 
gauge field theory in the 't Hooft large N number of colors limit [23]. To 
avoid unnecessary repetitions, we refer to the literature ED] • 

The reason why the above procedure works correctly for hermitian matrix 
models is the fact that the Green's function is guaranteed to be holomorphic 
in the whole complex plane except at most on one or more 1-dimensional 
intervals. 

The key difference which arises in the non-hermitian case (we denote 
the general non-hermitian matrix by X) is that eigenvalues of X are com- 
plex in general; in the large N limit they form two-dimensional domains 
in the complex plane, in contrary to one-dimensional cuts in the previous 
case. Therefore the power series expansion (jSJ) no longer captures the full 
information about the Green's function. In particular the eigenvalue dis- 
tribution is related to the non-analytic (non-holomorphic) behavior of the 
Green's function: 

-d- z G{z) = p(z) . (6) 

7T 

This phenomenon can be easily seen even in the the simplest non-hermitian 
ensemble — the Ginibre-Girko one [2311221, with non-hermitian matrices X, 
and measure 

P(X) = e ~ NTlXXt . (7) 

It is easy to verify that all moments vanish (trX™) = 0, for n > so the 
expansion (jSJ) gives the Green's function to be G(z) = 1/z. The true answer 
is, however, different. Only for \z\ > 1 one has indeed G(z) = 1/z. For 
\z\ < 1 the Green's function is nonholomorphic and equals G(z) = z. 

The above difficulty was first addressed in mathematical papers. Brown 
defined a measure for complex ensembles as 

= h (jm? + ww) logdet(x - x1n) { 



where 

det(X — Aljv) = exp 



^Trlog y^X-Al A )(.Yt-A:i , v , 



(9) 



is known in mathematics as Fuglede-Kadison determinant. For some recent 
results on Brown measure we refer to 
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Physicists have addressed the problem of measure, exploiting the analogy 
to two-dimensional electrostatics 0122 HZj • Let us define the "electrostatic 
potential" 

F = ^Tr\n[(zl N -X)(zl N -X^) + e 2 l N }. 

= ^lndet{(zl N -X)(zl N -X^) + e 2 l N ] (10) 



Then 



lim — — — - — - = lim — (Tr- 



6-0 dzdz e^oN\ (\zl N -X\ 2 +e 2 l N ) 2 

= ^(E^-aoWp(^) (ii) 

represents Gauss law, where z = x + iy. The last equality follows from the 
representation of the complex Dirac delta 



e 2 



7T 



^-A,) = lhn ^ ■ - (12 



*o (e 2 + \z - Aj 



In the spirit of the electrostatic analogy we can define the Green's function 
G(z,z), as an "electric field" 

G^J-WW *"- X \ 5 V (13) 

<9z iV^o\ (njv-Xt)^-^ +e 2 l Ar )/ v ; 

Then Gauss law leads (JEJ). 

The drawback of the above construction is that due to the quadratic 
structures in the denominator, form (|T3*jl is difficult to use in practical cal- 
culation. 

Instead of working ab initio with the object Ijlflj) . one can take consider 
the following generalization |2j. One defines the matrix- valued resolvent 
through 

Q = — /bTr ( z1n ~ X ielN 



with 
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N\ 2 I iel N zl N -X* 

n \""* \c d)) = \ g Tl q- x 

zl N - Xt 

(zl N -Xi)(zl N -X) + en N 

5 



B 

C 



—xe 

(zl N ~X)(zl N -X^) + en N 
—it 



(zl N - Xi)(zl N - X) + e 2 l 



N 



n _ zljy - X 

(zi N - x)(n N - xt) + en N { ] 

and where we introduced the 'block trace' defined as 

Then, by definition, the upper- right component Qu, hereafter denoted by 
Gx(z, z), is equal to the Green's function (fT3|) . 

The block approach has several advantages. First of all it is linear in 
the random matrices X allowing for a simple diagrammatic calculational 
procedure. Let us define 2N by 2N matrices 

Then the generalized Green's function is given formally by the same definition 
as the usual Green's function G, 

1 W 1 



n\ z e -n 

What is more important, also in this case the Green's function is completely 
determined by the knowledge of all matrix-valued moments 

<bTr 2 Z- X UZ- X U . . . Z- 1 ) . (19) 

This last observation allowed for a diagrammatic interpretation The 
Feynman rules were analogous to the hermitian ones, only now one has to 
keep track of the block structure of the matrices. 

Let us summarize the general properties 2\ of the matrix valued gener- 
alized Green's function. Each component of the matrix carries important 
information about the stochastic properties of the system. There are always 
two solutions for Qu(X), one holomorphic, another non-holomorphic. The 
second one leads, via Gauss law, to the eigenvalue distribution. The first one, 
the holomorphic one is not "spurious" , it represents the generating function 
for the real moments of the complex distribution j2H], but, by definition, 
cannot reproduce the non-holomorphic spectrum of the ensemble. 
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The shape of the "coastline" bordering the "sea" of complex eigenvalues is 
determined by the matching conditions for the two solutions, i.e. it is deter- 
mined by imposing on the non-holomorphic solution the condition Q12Q21 = 0. 
The off-diagonal elements have an interesting interpretation [29 . They rep- 
resent the correlator between the left and right eigenvectors (introduced in 

(^(L a \L a )(R a \R a )S(z-X a )\ = -^S 12 g 21 = -^C x {z,z) (20) 

On boundary of the domain of the eigenvalues, the above correlator vanishes. 

We would like to mention for completeness, that another, similar approach 
appears often in the literature, under the "hermitization method" name [SI 
IH E] . Basically it uses an alternative representation for the determinant in 
(|1U|). which is rewritten in a hermitian form 

det[(zl N -X)(zl N -X^)+e 2 l N ] = -det(^ 1n Yj zlN ~ X ) (21) 

Both versions lead to the similar results. 

An important feature of the explained above method of the generalized 
matrix-valued Green's functions is their link to Free Random Variables cal- 
culus [7J IB] , which we will address in the next section. 



3 Free Random Variables Approach 

In the hermitian RMT a method of free random variables (FRV), introduced 
by Voiculescu j7], is widely used ED E21 ESI- The power of the method 
lies in the fact, that in particular it allows to find the resolvent for the sum 
of two independent (free) 1 ensembles H% and H 2 , on the basis of knowledge 
of the resolvents for each separated ensemble only. In other words, oper- 
ation R introduced by Voiculescu linearizes the non-commutative, matrix 
convolution of hermitian random variables, alike the logarithm of character- 
istic function does it for sum of two identical, independent, one-dimensional 
random variables. 

In this paper we recall a version of Free Random Variables formulated 
by Zee jH] and known under the name Blue's functions. Blue's function is 
simply related to original R-function of Voiculescu 

B(z) = R{z) + - (22) 

^reeness is a non-commutative analog of statistical independence, introduced by 
Voiculescu 0. 
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and its usefulness comes from the fact, that it is simply the functional inverse 
of the resolvent 

B(G(z)) = G(B(z)) = z. (23) 

The algorithm of "adding" ensembles is therefore straightforward: consider 
two freely independent hermitian random matrices, H and H', and their 
(holomorphic) Green's functions, Gh{z) and Gh'{z). The addition algorithm 
goes as follows: First, one inverts functionally Gh(z) and Gh'(z) to get Blue's 
functions B H (z) and B H >(z). Second, using the addition law 

B H+H ,{z) = B H {z) + B H ,{z) - -. (24) 

z 

one obtains the Blue's function for the sum of ensembles. Third, inverting 
functionally this function one recovers needed Gh+h'(z). 

The generalization of the concept of the Blue's function for non-hermitian 
ensembles was first proposed in [Tj and confirmed by [HUH]. The generalized, 
matrix-valued Blue's function was defined as a functional inverse of matrix- 
valued Green's function, i.e. 

G(B(Z £ )) = Z e (25) 

for matrix Z (fTTj) . Such defined Blue's functions fulfilled the similar, but 
now matrix-valued addition law, i.e. 

B 1+2 (Z) = B X {Z) + B 2 {Z) - Z- 1 (26) 

with e safely put to zero. The explicit form of the matrix-valued Blue's 
function for an arbitrary ensembles was, however, far from obvious. In sev- 
eral applications using this method, an additional insight coming from e.g. 
diagrammatic interpretation was necessary. 

One of the purposes of this paper is to propose the rigorous method yield- 
ing the matrix-valued Blue's function solely on the basis of algebraic prop- 
erties, and abandoning any additional insight from e.g. diagrammatization 
methods. We will shown in the following chapters, that such construction is 
not only possible, but leads to strikingly simple algorithms for addition of a 
broad class of non-hermitian ensembles. 

4 Quaternion Green's functions 

The main advantage of working with the complex Green's function for her- 
mitian ensembles stems from the holomorphic properties of the resolvent. 
This function, being holomorphic everywhere except the intervals on real 
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axis, allows, by analytical calculation, to recover the spectral function by 
approaching the cuts from above and from below, i.e. by calculating the 
discontinuities along the cuts. 

It is tempting to propose the similar method in the case of complex spec- 
tra. A natural generalization is the algebra of quaternions. Such a spec- 
ulation appeared in the literature @], but the explicit realization of such 
construction was never completed and left as a challenge. In this chapter we 
propose the solution. We will now proceed in the following way: 

1. First, we introduce the notation. 

2. Second, we generalize the notion of a matrix-valued Green's function 
()18)) by exchanging Z e by a general quaternion Q thus exploiting the 
whole four-dimensional quaternion space and not only its part 'close' 
to the complex plane. A quaternion Blue's function is its functional 
inverse. 

3. Third, we compute quaternion Green's and Blue's functions for a her- 
mitian random matrix H. 

4. Fourth, we investigate how these functions behave during multiplying 
a general random matrix by a fixed complex number g. 

5. Fifth, we present an addition formula which enables to express a quater- 
nion Blue's function of a sum of two freely independent general random 
matrices by their respective quaternion Blue's functions. Then we com- 
pute a quaternion Blue's function for a non-hermitian X = H + iH', H 
and H' being two hermitian and freely independent random matrices. 

4.1 Notation 

We found it convenient to introduce the following notation: 
Writing down explicitly any square N x N matrix we add a subscript 'jy' to 
make the reader sure about the matrix' size. The blank places in a matrix 
mean that they are occupied by zeroes. We use standard Pauli matrices. If 
X is a general complex N x N matrix then 



where the superscript means 'duplication' since we duplicate X into X 
and X* . For quaternions, i. e. combinations of I2 and three Pauli matrices 
of the form 




(27) 



Q = x l 2 + ix ■ a, 
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we use a 2 x 2 matrix notation 



Q=(ai\. Mec (29) 

with a = Xq + 2x3, 6 = x\ + 2x2- If Q is a general quaternion, then 

Q u = f ^ ) , (30) 

where the superscript ' u ' means 'uniatal extension', since we extend Q by 
multiplying each of its elements by a unit matrix (N is a fixed number defined 
below) and also 

<? = «(' -Or " 73 "!" (31) 

where the superscript comes from imaginary units appearing in (j31j) . A 
quaternion of a special interest for us is defined when a complex variable z 
is being considered, 

*=(*,),■ <*> 

Generally, quaternion Green's (Blue's) functions we denote by calligraph 
alphabet to avoid confusion with complex standard Green's (Blue's) functions 
G{z) (B(z)). 

4.2 Definition of a Quaternion Green's Function 

Let us consider a random N x N matrix X, hermitian or not. We see that 
the key quantity for X is its matrix-valued Green's function (jlSjl : it has the 
form of a complex Green's function with a complex variable z substituted by 
a 2N x 2N matrix Z e u and X substituted by a 2N x 2N matrix X D . 

Let us consider a straightforward generalization of 1)180: from now on we 
call the quaternion Green's function the object 

Gx(Q) = ^<bTr 2 ^^>, (33) 
where Q is a quaternion 

a ib 



This defines as a quaternion function of a quaternion variable. 
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The only difference to the original meaning of this notion is that we 
exchanged Z e with a general quaternion Q thus promoting Qx to be the 
function of Q. Particularly, for Q — Z e we arrive at the original meaning pQ. 
The key idea of this work is that Qx has a few new interesting features as a 
function of Q which cannot be seen when regarding only the previous case 
of Q = Z t , i.e. assuming that e is infinitesimally small. 

The first interesting feature of Qx is that it enables to define the matrix 
Blue's function as the functional inverse of Qx'- 

Qx{B x {Q)) = B x (Qx(Q)) = Q, (35) 

for any Q. This is clearly the quaternion function of the quaternion variable. 
It is a natural generalization from the complex case (J2HJ)- Note that if Q is 
a quaternion, then Gx{Q) and Bx(Q) are quaternions too. It requires Qx as 
a function of Q everywhere and not only 'near' Z. It is particularly visible 
when studying the eigenvalues of Q. A quaternion Q (12~H|) has a very simple 
structure of eigenvalues, 

q x = q — x + i\x\, q 2 — q — x — i\x\. (36) 

The eigenvalues have following basic properties: 

• They are mutually conjugated. 

• They are different (there is no degeneracy) always except of a very 
special case of x = 0, which means that the only degenerate quaternion 
is the one of the form Q = Xq12- In other words, they are non-real 
always except the degenerate case. 

• \q\ 2 = %1 + x 2 = \a\ 2 + \b\ 2 = DetQ. Since the determinant is positive, 
it always fulfills the requirement of Fuglede-Kadison construction for 
Brown measure. 

Let us finish with speculative remarks. Introducing quaternions seems to 
be very natural from the algebraic point of view, however the reader will find 
by inspection, that for applications presented here not the whole quaternion 
space is necessary; it can be restricted to the subspace with x<i = 0. Since 
our calculations do not simplify considerably with regards to x 2 = case, 
we assume the general case x 2 ^ 0. Finally, we mention that the construc- 
tion presented here holds if instead of quaternions a general 2x2 matrix 
(biquaternion) is used. In this case the matrix-space is eight-dimensional, 
however some nice properties of quaternion eigenvalues are lost. 
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4.3 Quaternion Green's and Blue's functions for a Her- 
mitian Random Matrix 

Consider a known hermitian random N x N matrix H . We show how to 
calculate quaternion Green's and Blue's functions for hermitian H . The 
method is based upon a simple behavior of quaternion Green's functions un- 
der similarity transformations (jSH| and a knowledge of their form on diagonal 
matrices (J37j) . The same holds for Blue's functions, (jjDj) and (|3~S|). 

Let us first note, that a quaternion Green's function for H at the point of a 
diagonal matrix has a very simple form of a diagonal matrix with holomorphic 
Green's functions for H on the diagonal 

"wfM'-C^^l (37) 

Applying to both sides of (pT7j) and exchanging g, g with Gh(q), Gh(q) 
we get an identical theorem for a quaternion Blue's function for H: 

B 4\\H BM *«>),■ (38) 

To regain a general quaternion Q from its eigenvalues we use a similarity 
transformation S. The quaternion Green's function for H cooperates very 
well with similarity transformations of its argument; consider a similarity 
transformation S that is an internal operation in the quaternion space, then 

Qh^QS) = S- 1 g H (Q)S. (39) 

The similar formula holds for a quaternion Blue's function for H, 

BniS^QS) = S- 1 B H (Q)S. (40) 

Proofs are presented in appendix 1X1 

From these two properties we can also deduce that eigenvalues of Gh(Q) 
(resp. B H {Q)) are (G H (q),G H (q)) (resp. (B H (q), B H {q))). 

These two simple algebraic properties of quaternion Green's and Blue's 
functions, (pT7)). (J3~9"j) and (|3*5|). (jjOj) allow us to calculate them for any Q. Con- 
sider now any non-degenerate quaternion Q — \ \ - J with eigenvalues 



«6 a 

g, g and let S 1 be the similarity transformation diagonalizing Q, 



Q = *- l \ H q ) S. (41) 
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It is straightforward to apply expressions (}3*Tj) and (pl§j) to find a general form 
of QxiQ)- For a given Q we compute an explicit form of S and S~ l , then a 
simple three-matrix multiplication leads to the result. The result is derived 
in appendix iBl 

Q H {Q)=lH{q,q)l2-i H {q,q)Q\ (42) 

where jh and j' H are two scalar functions depending only on Q's eigenvalues 
and are given by 

iu(qi) = «G»fa)-gC»(fl (43) 
, = G "M~ G ^ . (44) 

q-q 

The degenerate case of Q = X0I2 is trivial since the degenerate quaternion is 
already diagonal, so that from (JHZj) 

fe(x l 2 ) =Gh(xo)1 2 . (45) 

This is the solution; we have expressed Qh by Gh, i-e. the quaternion gen- 
eralization using the known, holomorphic resolvent. 

Let us note, that due to (jHHjl and PH|) the identical construction can be 
made for Bh', hence the analogous expression for non-degenerate quaternion 

Q, 

B H {Q) = Ph(Q, Q)h - P H (q, q)Q\ (46) 
where (3h and (3' H are similar to -yn and j' H but with B H in the place of Gh, 

PhM = -fl>»(fl (47) 

#,(?,?) = (48) 
q-q 

If a degeneracy is present, a formula analogous to holds. Again, we have 
expressed quaternion by standard analytic Bh- 

4.4 Quaternion Green's and Blue's functions for a Gen- 
eral Random Matrix Multiplied by a Fixed Com- 
plex Number 

We have already computed quaternion Green's and Blue's functions for a 
hermitian random matrix H. Now there is time for the passage to a non- 
hermitian random matrix. The following problem is the introductory one 
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before doing so: consider a general random matrix X and a fixed complex 
number g. We want to express the quaternion Green's and Blue's functions 
for gX through the ones for X. 

If g = then trivially Qo(Q) = ^, so let us further assume that g ^ 0. In 
appendix O we present a simple derivation of 

» = w ( 1/9 i/J 2 e) ( V9 i/s),- (49) 



This would be incorrect unless the argument were just a quaternion Q; 
1/9 . . 

a — > a/gf, 6 — >■ fe/g (50) 



in fact I ^'^ ,_ I Q is a quaternion with 



(see (12911 ). In terms of variables (JZHj), this replacement is non-trivial. 

For a special case of real g, Eq. (flTij) simplifies, 

6U(Q) = -Sx(-Ql geR\{0}. (51) 

Let us note that this formula is identical to the one well known for a complex 
Green's function and complex g, 

G gX {z) = -G x {-z). (52) 
9 9 

Similar behavior appears at the level of Blue's functions. Again, if g = 0, 
then trivially Bq(Q) = 4, so let us further assume that g ^ 0. From the 
definition and (|4*H|) we immediately get 

«W« = ( 9 j), B *«( 9 s ),>- P3) 

The argument here is a quaternion with 

a — > ag, b — >■ 6^. (54) 

For a special case of real g Eq. (}53|) simplifies, 

B gX {Q) = gB x (gQ), geR\{0}. (55) 

Let us note that this formula is identical to the one for a complex Blue's 
function and complex g, 

B gX (z)=gB x (gz). (56) 
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4.5 A Quaternion Blue's Function for a Non— Hermit ian 
Random Matrix with Hermitian and Anti— Hermitian 
Parts Freely Independent 

In this subsection we apply a general formalism of the previous part to a 
particular problem: consider a non hermitian random N x N matrix X of 
the form 

X = H + iH', (57) 

where the (hermitian) random matrices H and H' are freely independent. 
The advantage of quaternion Blue's functions stems from the fact, that, not- 
surprisingly, they obey the addition law (the proof parallels the proof of an 
"addition law" for hermitized matrices in [3], so we are not repeating it here) 

B X (Q) = B H (Q) + BiH'{Q) - ^- (58) 

We can now easily compute the quaternion Blue's function for X assuming 
the knowledge of quaternion Blue's functions for H and H', since they are 
expressed using (jUj) , i.e. their complex Blue functions. 

Using pUj) and ([53)1 . we may rewrite (|58|) in an explicit form, 



B x (Q) = (3 H (q,q)l 2 + (3 H/ (q 1 ,q 



i TIN ' 1 



- W H (q, q) + /W, Q 1 ) + v^qW' ( 59 ) 

where q 1 , q 1 are eigenvalues of Q l , see fl3*T|) . 

This explicit form of the addition law for any ensemble of the form H+iH 
(where H and H are free) is one of the main results of this paper. 

Few comments are helpful. 

• The : 's terms came from the multiplication formula ()53|) and the non- 
real number i before H . Note that the presence of % in (|H7j) is essential 
for the non-hermiticity and it is % that brings Q l , so the existence of 
the eigenvalues of Q l on top of eigenvalues of Q is exactly a sign of 
non-hermiticity. 

• The operation 1 changes the trace of Q and does not change its deter- 
minant, 

TrQ 1 = i(a - a) = -2x 3 ^ 2x = TrQ, DetQ 1 = DetQ, (60) 
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see (J2H), so that TrQ 1 cannot be expressed through Q's invariants. 
Since in general 



q = x + i\x\ = x + iyxj + x\ + x|, (61) 

then 

q 1 = -x 3 + %\jx\ + x\ + x%, (62) 

and it is impossible to express q l only by q. Note that the eigenvalues 
are related since 

\q\ = \q l \. (63) 

• The algebraic properties (fHTj) . (jHEJ) of quaternion Green's and Blue's 
functions for hermitian random matrices remain valid for general ran- 
dom matrices, 

Bx{q) \ ( B x {q) \ 

Bxt(q)J 2 V B x (q) ) 2 > 

(64) 

and identically for Q X - (We use general formulae, B x ^(q) = B x (q) and 
identically for G.) It is important to note, that due to the presence of 

B x {S- l QS) ^ S- 1 B X (Q)S (65) 

and similar relation holds for Q x . Indeed, this important non-trivial 
behavior under the similarity transformation is precisely a footprint of 
non-hermiticity of X. The immediate conclusion is that (B x (q) and 
B x (q)) (resp. (G x (q) and G x (q))) are not the eigenvalues of B X {Q) 
(resp. Qx{Q))- 

5 Holomorphic and Non— Holomorphic Green's 
Functions for a Non— Hermitian Random 
Matrix with Hermitian and Anti— Hermitian 
Parts Freely Independent 

In this section, we invert functionally in a fixed point of Q = Z which gives 
Q X (Z) and the (complex) non-holomorphic Green's function for X in par- 
ticular. The result is given by a surprisingly simple and operationally con- 
venient formula. We have already computed the quaternion Blue's function 
for X = H + iH . We should now invert it functionally to get the quaternion 
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Green function for X, however this is in fact unnecessary since in this paper 
we are interested only in Q X {Z). Hence we need only to solve the equation 

*<U "WO.- (66) 

[A iB \ 

where we denote Qx(Z) = I . ^ t ] . We use it to extract A = Gx(z, z) 

\ %B A J 2 

and — \B\ 2 = Cx(z, z). 

One may ask why it is allowed to put here the regulator e = 0. Note 
first that QxiQ) is discontinuous in Q = Z, i. e. Qx{Z) ^ Gx(Z e )- On 
the other hand, we will show that there are always two solutions of (|66|. 
the first one gives Qx(Z) and the second one, Qx(Z e ). Finally, the answer 
is that functional inverting of a function in her discontinuity point gives two 
solutions, the first one is a value of the function in this point and the second 
one is a limit value of the function when approaching the discontinuity point. 
In other word, the limit e — > is precisely encoded in the functional inverting. 

Rewrite using (JoTJj) . 

k + ik' — lA HB 

UB k-ik'-lA )J~ V z 

where for short 



(67) 



k = p H (Gx(Z)), (68) 
k' = (3 H >(Gx(Z) 1 ), (69) 

i = (3> H (g x (z)) + (3> Hf (gx(z?) + 1 (70) 

Detg x {Z) 

Looking at the two off-diagonal equations we see that this equation has 
always two solutions, one with B = and one with B ^ 0. 

Let us infer the first case. For B = 0, Qx(Z) is diagonal, hence from (jHljl 

'M A AH™ ™\- (71) 

which leads to 

A = G x (z). (72) 

We may call this solution the holomorphic one since it is expressed by a 
holomorphic Green's function for X. 

The second case is a general one and we call it the non-holomorphic one 
since then A is a non-holomorphic Green's function for X. 

To sum it up, there always exist two solutions of A 6' 6)) . the first one with 
B = 0, which reads A = Gx(z) and which is called the holomorphic solution 
and the second one with B ^ which is called the non-holomorphic solution. 
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5.1 The Non— Holomorphic Solution 

Since in (}68j) . (}69|) . (fTU|) functions /3 and /3' are expressed via eigenvalues of 
Qx(Z) and Gx(Z) 1 , let us denote these eigenvalues by g and gr 1 , and their 
conjugate partners respectively. From subsection (|4.5jl we know that they 
are not equal to Gx(z), Gx(z). 

In the non-holomorphic case the off-diagonal equations of ()66|) are 



/ = 0, (73) 
because it is possible to divide both sides by iB or %B\ explicitly 
B H {g) - B H (g) + B H ,{g l ) - B H ,@) + J_ = Q 

9-9 g l -g l \9? 

Two diagonal equations of (j66|l 

k + ik' — ID = z, 
k — ik' — I A = z 



gives due to f73^l 



k = x, (75) 
k' = y, (76) 



where z = x + iy. Explicitly 

gB H (g) -gB H (g) 
9-9 

g l B H ,(g l )-jB HI (J) 



x, (77) 
V- (78) 



The form of these equations suggests introducing two new (complex in 
general) variables, m and m', thus changing two equations (|77jl . (|7Hj) into 
four equations 

gB H {g) = xg + m, 

gB H {.g) = xg + m, 
jB H ,tf) = yg^ + m', 
g l B H ,(g l ) = yg l + m 

Conjugating the first and third ones, comparing with the remaining two and 
exploiting B H (z) = B H {z) true for hermitian matrices, we get 

m,m eK (79) 
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and two equations 



Tfl 

B H {g) = x+-, (80) 
9 

777 ' 

B H <{g l ) = y + —. (81) 
9 

Substituting (jEDJ), (jHTjl into the off-diagonal equation (|7Ijl gives a simple 
result, 

m + m' = l. (82) 

To bring the solution to the end we have to express A and B through g 
and g 1 ; using general expressions (JfiTf and (J62|) we get 

A = 9fy - (83) 

and 

- \B\ 2 = l -m 9 f + (^) 2 ) - ((%) 2 + (% 1 ) 2 )]- (84) 

Finally, we can write down the main result of this section: the algorithm 
of calculating a non-holomorphic Green function for a non-hermitian random 
matrix X with hermitian and anti-hermitian parts freely independent. 

Let us start with given non-hermitian random matrix X = H + iH', with 
H and H' freely independent hermitian ensembles. The holomorphic Green's 
or Blue's functions for H and H' are explicitly known or given implicitly by 
some equations. The algorithm goes as follows: 

1. Write down two equations 

Tfl 

B H (g) = x + - (85) 
9 

1 YY) 

B H >{g l ) = y + ^> (86) 

with three unknown quantities, g, g l and m G M. Find from them g + g, 
gg-i g 1 + g l an d g l g l expressed via m. 

2. Compute m from the third equation, 

99 = 9 l ¥- (87) 

3. Put m into g + g, g l + g l and \g\ 2 expressed through m (step 1). 
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These three steps yield the non-holomorphic Green's function and corre- 



lator between left and right eigenvectors for X: 

Gx(x,y) = \[{g + g)-i{g l + J 1 )], (88) 

C x (x,y) = \(g 2 + f + (g 1 ) 2 + g?) = (89) 

= \[{g + g) 2 + {g l + 7f\-\g\ 2 - (90) 



The general algorithm contains also a simple method of deriving an equation 
(in coordinates (x, y)) of a borderline of X's eigenvalues' domains. The bor- 
derline is exactly the place where holomorphic and non-holomorphic Green 
functions meet together. And it is the limit B — > that carries the non- 
holomorphic solution towards its boundary with the holomorphic one. Hence 
the borderline 's equation is simply 

5 = 0, (91) 

where B is an appropriate element of the non-holomorphic solution. This 
is equivalent to the condition of vanishing of the correlator between left and 
right eigenvectors, 

C x (x,y)=0, (92) 

where we compute Cx {x, y) for the non-holomorphic solution. 

Due to (J90|) and ijBTjl . the general equation for the curve defining the 
support of the eigenvalue domains is given by 

^ + 9) 2 + {i + 7) 2 = ^99, (93) 

where g + g, g l + g 1 and gg are expressed via m in step 1 and explicitly in 
step 3, so that this equation binds x and y to form the borderline. 

Before closing this section, few comments are useful. First, inferring the 
equations (|85|l . (JHHJ) we look for a solution that comes with its conjugate 
partner in pair. (E. g., if g satisfies the equation, g does it as well.) So that 
two possibilities are allowed, either this is the pair of two different (and 
hence non-real) numbers, or they merge into a one real number. The second 
possibility means actually, that Qx{Z) is proportional (with real coefficient) 
to a diagonal matrix, i. e. (a) in case of g, eq. (|55J). to 1 2 (b) in case of 
g , eq. (J86|) . to 203. However, this is impossible since we are looking for 
a non-holomorphic solution with non-diagonal Qx{Z). To sum up, we are 
looking for a solution that is non-real and comes in pair with its conjugate 
partner. 
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Second, it is not actually our task to find all g's in this step but only their 
sums, products etc. , e.g. we may use the Viete rules. 

Third, note that it is possible to use the holomorphic Green's functions 
for H or H' rather than Blue's ones in (|85j) . (jHEJ), by inverting these equations 
functionally. In some cases this approach may be considerably simpler. 

6 The Version of the Algorithm in the Case 
of H from GUE 

A common case is when one of the ensembles belongs to Gaussian class, 
but for a second a holomorphic Blue's function is given by a complicated 
equation. In this case one may reformulate slightly the algorithm. 
Here let us deal with the problem of Gaussian H with 

B H (s) = s + - (94) 
s 

and very complicated B H /(s). If so, equations (JS3J), (|SS]l are very complicated 
and another way is preferred. 

The quaternion Blue's function for H is obtained easily from (|4fi)l and 
reads 

B H (Q) = Q + ± 

hence the equation (|66p. after short manipulation, becomes 

/ iA B \ ( i(A-z) -B 
Bh '\ -B -iA )J-{ B i(z-A) 

The key observation is to invert the last equation functionally. 

Let h, h be the eigenvalues of ( ^„ ^ .,_ j , i. e. from the 

^ B i(z-A) J 2 ' 

general formula 

h = y - + i v /{$lA-x) 2 + \B\ 2 ), 

where we recall z = x + iy. 
Hence from 



-B * a )=^ h ~ h ^-^ h ^( i( *-B A) i(A B -z) ) < 95 > 
which is the equation to solve. 
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The off-diagonal equations of ()95|) are thus simply 

' fu m G H ,{h) - G H '{h) 
j H ,(h, h) = — j = -1. (96) 

The diagonal ones 

iA = 7h'(/i, h) - j' H ,(h, h)i(z - A), 

become due to (j96|) 

iA = j H '{h, h)+i{z-A), 

or equivalent ly 

- hG H ,(h)-hG H ,(h) 

nf H '{h,h) = — = = -y, (97) 

h — h 

A + A = x. (98) 

Note differences between @J7|), JSHJ) and 1(77)1 . (I7H1) . We solve the second 
diagonal equation, (j9*8~|) . by introducing a new unknown quantity n, 

A = ^+n. (99) 

The off-diagonal equation (jHEJ) is solved by introducing a new unknown vari- 
able m, 

G H ,(h) = -h + m, meR. (100) 
Inserting this into the first diagonal equation (jHTJ we get easily m, 

m = h + h-y; (101) 

but h + h = 2y — 2n, so that 

y — m 

m = y — 2n — > n = 

and hence 

A=^p. (102) 

Finally, it is now possible to write down the algorithm of calculating a non- 
holomorphic Green's function for a non-hermitian random matrix X with 
hermitian and anti-hermitian parts freely independent and with hermitian 
part being a GUE random matrix. 

We consider a non-hermitian random matrix X = H+iH', with H and H' 
freely independent and Bh(s) = s + ^, is given. The complex holomorphic 
Green's or Blue functions for H' is known (or more often) given by some 
(complicated) equation. The algorithm reads: 
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1. Write down an equation 



G H ,(h) = -h + m, (103) 

with two unknown quantities, h and m G R. Find from it h + h and 
hh expressed via m. 

2. Compute m from the second equation, 

h + h = y + m. (104) 

3. (Necessary only for Cx and the borderline's equation.) Put m into the 
expression for hh. 

The results of these three steps are given explicitly by 

Gx(x,y) = Z -^- , (105) 
C x (x,y) = ~[x 2 + (y + m) 2 } - hh. (106) 

The borderline's equation is 

Ahh = x 2 + (y + m) 2 , (107) 
where m is known from step 2 and hh from step 3. 



6.1 Examples 

To demonstrate the usefulness of the advocated approach, we recalculate 
three classical examples well known in the literature. We are not describing 
these models in detail, neither their physical context. We start with known 
holomorphic Green's (or Blue's) functions for hermitian and anti-hermitian 
parts, referring for further details to the original papers. 

• The Girko-Ginibre Model |2H E2J - 

Gaussian complex ensemble, known also as a Girko-Ginibre ensemble, corre- 
sponds to the case when X = H + iH , where both H and H are Gaussian 
Unitary Ensembles. Since the holomorphic Green's function in this case 
equals simply to G(z) = (z — \/ z 2 — 4)/2, with imaginary part yielding a 
seminal Wigner semicircle [34j, the functional inverses of the resolvents are 
simply 

B H (s) = rs + -, B H >(s) = r's + -, (108) 
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where r, r' are real and positive constants (arbitrary variances of the Gaussian 
distributions). We apply the general algorithm: 
Step 1: Equations (JEEJ), (J86j) are quadratic 

rg 2 — xg + 1 — m = 0, r'(g 1 ) 2 — y(g l ) + m = 0, 

hence from Viete rules 

x 2 1 — m 

9 + 9 = ~, \g\ = , 

I , — 2/ | 1,2 m 



Step 2: The equation for m is 



1 — m m r' 

m 



fy qyl qy j iyf 

Step 3: We get finally 



Gx(x,y) = \{- r -i^), (109) 
Cx(x,y) = \(^ + ^)--±- i (HO) 



or in the special case of r = r' = 1, 



Cx(x,y) = ~z, (111) 
C*(*,l/) = ^ 2 + 2/ 2 )-i (112) 



Borderline is given by (|11U|) 



x 2 y 2 4 



(113) 



which is the ellipse with semi-axes -^grpa an d y/X+r^ ' ^ n ^ e case °^ 
r = r' = 1 this is a circle with radius \/2. Applying Gauss law to non- 
holomphic solution yields spectral distribution (here a constant). 

•Model of the Chaotic Resonance Scattering [TH] 
Non-hermitian ensemble X = H + iH is defined by H belonging to GUE (or 
GOE) and if is a Wishart ensemble. Corresponding Blue's functions were 
e.g. calculated in [Tj l2Hj. 

I CT I 

B H ( s ) = s + -, B w (s) = -—— + -, (114) 
s 1 + cs s 
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where r, c are real constants. We apply the general algorithm. 
Step 1: 

g 2 — xg + 1 — m = 0, yc(g 1 ) 2 + (cr — cm + y)g l — m = 0, 
hence from Viete rules 

g + g = x, \g\ 2 = l-m, 

J|2 



T - cm — cr — y 
9 l + 9 l = 



yc 



Step 2: The equation for m is 

m 



1 — m = 



\9 



m 



m 
yc' 



yc 



yc 



yc-1 



Step 3 gives 



g + g = x, g l + g l 



r 1 



yc-1 y c 



1 - yc 



thence 



G x (x,y) = -[x-i( — ^— ----- 
2 yc — 1 y c 

r 1 



)], 



C^(x,y) = \[x 2 + {^-----) 2 ]--^. 

4 yc— 1 y c 1 — yc 



The borderline's equation is 



x 2 + 



r--) 2 = 



yc — 1 y c 1 — yc 



(115) 
(116) 

(117) 



•The Complex Pastur Model 

By complex Pastur we mean a class X = H + iH\ where H is deterministic 
and H' belongs to GUE. We choose for simplicity the deterministic part with 
only two distinct opposite eigenvalues, so 



G N (s) 



1 1 

+ 



S — [l S + jJ, 



(118) 



where \x is a real constant. The general algorithm is being applied. 
Step 1: 



* + 7 

9 (X+S)2_ 



(x 2 - /i 2 )y 2 + i(2m - l)x + m(m - 1) = 0, 
25 



{g l Y -yg 1 + m = o, 

hence from Viete rules 

x(l — 2m) l9 m(m — 1) 
9 + 9 = — 5 J 2 ' # =— 1 2^" 

g l + 9 I = V, l^ 1 ! 2 = m. 
Step 2: The equation for m reads 

m(m — 1) 



x 2 — fi 2 



m — ► m = x 2 — /i 2 + 1 or m — 0. 



Step 3 gives 

+ = -2a; — 2 x , , g l + g l = y, \g\ 2 = x 2 - /j 2 + l 

or 

# + # = ^ 2> 9 1 +9 I = V, l^| 2 = - 

In the second case the first and third equation contradict each other, so that 
finally 

Gx(x,y) = -{-*- 2(x2 X _^y (H9) 
Cx(x,y) = \[{2x + * — f + y 2 ] - (x 2 - ^ 2 + 1). (120) 

The borderline's equation reads 

y 2 = 4(* 2 - /i 2 + 1) - (2x + 2{x2 X _^ 2) ) 2 (121) 

or after some manipulations reducing it to the form calculated first by Stephanov 



y 2 = -r-~ ■ - -r- ■ -r- v- r- , _ (m) 



-4//V + x 2 (8/i 4 - V - 1) + 4/i 4 (l - /2 



2\ 



(x 2 — jJ?) 2 



2 Original calculation involved chiral GUE, but in the leading N one-point Green's func- 
tion does not feel the difference in the bulk of the spectra. Note, however, that two-point 
Green's functions are different for chiral and non-chiral complex Pastur ensembles 35 . 
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7 Summary 



In this paper we introduced the concept of quaternion Green's function for 
non-hermitian ensembles with complex spectra, borrowing from the analogy 
of complex Green's function for hermitian ensembles with real spectra. The 
resulting matrix resolvent encodes important spectral properties. First, there 
are always two solutions for the diagonal elements (e.g. Qn). The nonholo- 
morphic one yields, via Gauss law, the spectral density. The second solution 
is holomorphic. Both solutions match along the spectral curve, defining the 
support of the eigenvalues. The product of off-diagonal elements gives 
the correlator between the left and right eigenvectors, hence the shape of the 
spectral curve can be also inferred by equating the value of this correlator 
to zero. The functional inverses of the quaternion Green's function, named 
quaternion Blue's function, obey for statistically independent ensembles an 
addition law, following precisely similar laws in Free Random Variables cal- 
culus for the case of hermitian ensembles. All these results confirm and 
generalize previous observations made recently in the literature. In partic- 
ular, construction [TJ |2] can be viewed as a special case of the quaternion. 
The similar construction proposed in [HI EJ E] can be viewed as a variant 
of the presented here, with quaternion replaced by a two-by-two hermitian 
matrix. The advantage of the approach presented here lies in the fact, that 
the defined procedure is purely algebraic, and does not involve any additional 
implicit insights from e.g. diagrammatic approaches, used in the other works. 
It is also simpler at the operational level then the hermitization procedure. 
Last but not least, we dare to say that the quaternion extension in the case 
of complex spectra is an esthetic construction. 

We proposed the general, simple, algebraic algorithm of finding the matrix- 
valued Green's function for an arbitrary non-hermitian random matrix mod- 
els with hermitian and anti-hermitian parts freely independent. This reduced 
the problem to solving some given explicitly (usually polynomial) equations. 
The resulting algorithm is surprisingly simple, as we hope the reader does in- 
fer by comparing three classical examples recalculated here using our method 
to original derivations. 

We concentrated here on mathematical aspects and the clarity of the pre- 
sentation, postponing announcing new results for particular physical models 
to the approaching papers. 

We considered here only random non-hermitian ensembles with complex 
entries, not discussing the real or quaternion ensembles. In the case of one- 
point function and in the leading large N expansion, the Green's functions 
for similar real or quaternion cases are expected to be similar modulo triv- 
ial rescalings. However, challenging problem is the subleading behavior, or 
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two-point Green's functions (wide-correlators) for nonhermitian ensembles. 
We expect that an extension of the quaternion technique may considerably 
simplify the calculations of wide-correlators for two- and more-point Green's 
functions for non-hermitian ensembles. These ideas will be addressed next. 
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A Behavior of Quaternion Green's and Blue's 
Functions for a Hermitian Random Matrix 
under Similarity Transformations 

In this appendix we present proofs of theorems (|3^|) and (fIUj). 

A.l Proof of ( 1891) 

Since H is hermitian, 

H°=( H ) . (123) 

V / 27V 

The diagonal structure and the equality of the diagonal elements implies that 
H D commutes with any matrix of the form 5* u , 

[H D ,S v }=0; (124) 

it is easy to check. Hence 

g-l,V H T> S U = R D ( 125 ) 

Now let us choose S to be a similarity transformation that is internal in 
the quaternion space, i. e. if Q is a quaternion, then S~ 1 QS too. There exist 
such transformations as is shown in appendix iBl 

Now from the definition (}3*3*|) and this simple property of H D , ()125j) . we 

have 

GniS-'QS) = l(bTr 2 [(5- 1 g5) u -iJ D ]- 1 ) = 

= hbTr 2 [S- 1 ' v (Q v -H D )S l] }~ 1 ) = ... 
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The S v s can be get out of the block-trace (and the average of course) because 
of the simple relation for a block-trace, 

bTr 2 (FM u ) = (bTr 2 F)M (126) 

true for any 2N x 2N matrix Y and 2x2 matrix M. Hence 

. . . = S-^bTr^Q 11 - H^S = S- 1 g H (Q)S, 
what was to prove. 

A. 2 Proof of (HOD 

From ()39)1 and from the definition (}35|) we have 

g H {s- l QS) = s- 1 g H (Q)s 

s-'qs = B H (S- 1 g H (Q)S) 
S- 1 B H (Q)S = B H {S- l g H (B H {Q))S) = b h {s- x qs), 
what was to prove. 

B Quaternion Green's and Blue's Functions 
for a Hermitian Random Matrix 

In this appendix we present proofs of theorems and and also some 
additional information. 



B.l Proof of (H2D and 



The proof comes in three parts, first for a given Q we compute the similarity 
transformation S which diagonalizes Q, second we invert it to get S 1-1 and 
third we use these explicit forms to apply (|57j) and (|39"j) to get the final result. 

Let us assume, that Q — ( °1 ) is non-degenerate (eigenvalues are 

\ ib a J 

q 7^ q) since we know the degenerate case be trivial, Q = Xq12- Let us assume 
also that Q is not diagonal, which means b ^ 0, since the opposite case does 
not need any diagonalizing transformation. 

First let us find an explicit expression for S assuming the knowledge of a 
and b. A short calculation shows that 



S 



a—q 

Si si-- 



a—q 



b 
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for any complex s±, $2- This is a general similarity transformation that is 
internal in the quaternion space. Let us choose for simplicity e. g. S\ = ib, 
S2 = ib, which gives 

S=( ib q 7r a ) . (127) 
\ q - a ib J 2 

(The 21 element q — a can be exchanged by a — q due to a + a = TrQ = q + q.) 
A determinant reads after a simple calculation 

DetS = (q-a)(q-q); (128) 

it is always non-zero (S is invertible) in the considered situation. 
The expression for S^ 1 will obviously be needed; easily 

S- 1 = — ( V ^ ) " ( 129 ) 
a — a \ 1 — / 



q-q 

Finally 



q—a 



2 



Qh{Q)=Qh{S- 1 [ q s ) S) 



Q 



2 



1 (A - 1 W G H (q) 

ib 



<!-<A I AAV 

z6 9 — a 
g — a zfe 



2 



Noticing also — \b\ 2 = (a — q)(a — q) and multiplying these three matrices we 
immediately recover the final result 



_ qG H (q) ~ qG H (q) 1 G H (q) ~ G H (q) nj 

... — = 1 2 = . 

q-q q-q 

The identical construction leads to (l46|) . 

B.2 Some Expressions Helpful to Cross— Check 
and 



Let us note the following easy to check relations, 



MB H «,),B Hm = |||, (130) 
i H (B Mh B H m = _L__ (131) 
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and similarly in the opposite way. 
Moreover, all 7s and /3s are real. 

Using these expressions we can e. g. immediately show the cross-check 
relation Gb(&h(Q)) = &h(Gh(Q)) = Q which ensures us that the derivation 
is correct. 

C Quaternion Green's and Blue's Functions 
for a General Random Matrix Multiplied 
by a Fixed Complex Number 

From the definition ()33j) and for complex g 7^ 0: 



Q 9 x{Q) = ^(bTr^QU - (gXr)- 1 ) = 
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